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CROSS-REFERENCE TO RELATED APPLICATION(S) 

[1001] This application claims the benefit under 35 U.S.C. § 1 19(e) of application 
number 60/432,024, filed December 10, 2002, which application is incorporated 
herein by reference. 

BACKGROUND 
Field of the Invention 

[1002] The invention relates to contact sensitive devices. 
Description of the Related Art 

[1003] Visual displays often include some form of touch sensitive screen. This is 
becoming more common with the emergence of the next generation of portable 
multimedia devices such as palm top computers. The most established technology 
using waves to detect contact is Surface Acoustic Wave (SAW), which generates high 
frequency waves on the surface of a glass screen, and their attenuation by the contact 
of a finger is used to detect the touch location. This technique is "time-of-flight", 
where the time for the disturbance to reach one or more sensors is used to detect the 
location. Such an approach is possible when the medium behaves in a non-dispersive 
manner i.e. the velocity of the waves does not vary significantly over the frequency 
range of interest. 



SUMMARY 

[1004] According to one embodiment of the invention, there is provided a contact 
sensitive device comprising a member capable of supporting bending waves, a 
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plurality (e.g. three or more) sensors mounted on the member for measuring bending 
wave vibration in the member, whereby each sensor determines a measured bending 
wave signal . A processor calculates a location of a contact on the member from the 
measured bending wave signals, in that the processor calculates a phase angle for each 
measured bending wave signal, then calculates a phase difference between the phase 
angles of least two pairs of sensors from which the location of the contact is 
determined. 

[1005] According to another embodiment of the invention, there is provided a 
method of determining information relating to a contact on a contact sensitive device 
that includes providing a member capable of supporting bending waves and a plurality 
of sensors (e.g., three or more) mounted on the member for measuring bending wave 
vibration in the member, applying a contact to the member at a location, using each 
sensor to determine a measured bending wave signal and calculating the location of a 
contact from the measured bending wave signal by calculating a phase angle for each 
measured bending wave signal, calculating a phase difference between the phase 
angles of at least two pairs of sensors and determining the location of the contact from 
the at least two calculated phase differences. 

[1006] The following features may be applied to both the device and the method 
with the processor being adapted to provide many of the calculations or processing 
steps of the method. 

[1007] Reflected waves may be suppressed by placing an absorber in contact with 
the edges of the member. The mechanical impedance of the absorber and member 
may be selected so as to significantly reduce or minimize reflections of bending 
waves from the edges of the member. In particular, the impedances may be selected 
such that bending wave energy is strongly absorbed in a frequency band around a 
chosen frequency co 0 . The impedance of the absorber may be selected to be both 
resistive and compliant. The impedances may be selected to satisfy the following 
equation: 

Z T =-iZ B (co 0 ) 

where Z T is the termination impedance of the absorber, Z B is the mechanical 
impedance of the edge of the member, and i is the imaginary number = V-T . 
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[1008] The absorber may be made from foamed plastics which may have open or 
closed cells and may be polyurethane or polyvinylchloride. For example the foam 
may be a soft PVC predominantly closed cell foam such as MIERS™ or a medium to 
high density, open cell polyurethane foam. The absorber may extend substantially 
around the periphery of the member. The absorber may act as a mounting which 
supports the member in a frame or to another surface. 

[1009] The member may comprise a raised pattern on its surface whereby a 
contact drawn across the surface provides a variable force to the member to generate 
bending waves in the member. The pattern may be periodic, or quasi-periodic with a 
statistically well-defined spatial distribution of undulations. The pattern may be 
random whereby a contact travelling over the surface of the member generates a 
random bending wave signal. The random relief pattern may be an anti-reflective 
coating, an anti-glare surface finish, or an etched finish such as those which are found 
on many known transparent panels placed in front of electronic displays. 

[1010] The member may be a liquid crystal display screen comprising liquid 
crystals utilized to excite or sense bending wave vibration in the member. 

[1011] Each measured bending wave signal may be processed by a band-pass 
filter with a pass-band centered at the chosen frequency coo and having a bandwidth of 
Aco. The bandwidth Aco of the filter is preferably chosen to address the Doppler effect 
whereby a bending wave arrives at a point with a different frequency from its original 
frequency. Accordingly, the bandwidth preferably obeys the relationship: 
Aco » 2£(<y 0 )v niax 

where v max is the maximum lateral velocity of the contact across the surface, e.g. if the 
contact is provided by a stylus, v max is the maximum velocity that a user is capable of 
moving the stylus. 

[1012] The phase of each filtered signal may be measured by comparison with a 
reference signal. The reference signal may have a frequency coo* The measured phase 
is the average phase difference between the input and reference signals, optimally 
measured over the interval 27i/Aco. Alternatively, the reference signal may be derived 



Final patent application 08l-003.doc 
Client Reference: NXT 222 



-3 - 



Attorney Docket No.: 081-0003 



from a filtered signal from a second sensor in which case, the measured phase is the 
phase difference between two input signals. 

[1013] The phase differences may be calculated at intervals of 27r/Aco which may 
be intervals of less than 10ms. The reference and input signals may be fed to a phase 
detector. Output from the phase detectors may be fed through low-pass filters having 
frequency cut-offs of approximately Aca/2, then through digitisers and finally through 
a processor to calculate the phase angle ft 

[1014] The instantaneous phases, and 9 m (t), of two measured bending wave 
signals may satisfy the phase difference equation: 
*0i m = 0/ " 0 m = k(co 0 )Ax iip + 2m lm 
where Ax lm = jc, - x m (x m and jc/ being the distance from the contact location to each 

sensor labelled m and 1 respectively), and £(co) is the wavevector. This equation may 
be satisfied if the path length difference between two sensors is less than the 
coherence length of the bandpass filter, which is defined as 

__ 2tvco q 

The coherence condition is therefore \Ax lm \ «x c . If the coherence condition is not 
satisfied, the above phase equation may not be satisfied. 

[1015] Thus, values of n im and the phase angle difference are required to 
determine the location of the contact. The shape of the member may be selected to 
constrain the magnitude of Ax im to values less than half of one wavelength, ie. 
|Ax /m | < 7r/k(a> 0 ). In this case, where all possible values of Ax im satisfy the condition 

|Ax /m | < x/k(a> 0 ) 9 there is only one value of n/ m which is the integer n im satisfying 

|A0 /m - 27m lm \ < n . Alternatively, n may be estimated or inferred in some way. 

[1016] Another class of foams that has been found to be suitable are acrylic 
closed cell foams. These may have a high degree of damping and relatively high 
stiffness. Such properties are particularly suited to edge termination of stiff, heavy 
materials such as glass. Examples include 3M serial numbers 4956, 4910, 4950, and 
4655. 
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[1017] Each phase angle difference in combination with the range of possible 
values of the integer n im may be used to generate a series of path length differences 
thereby defining a series of discrete hyperbolic curves on the surface of the member, 
denoting possible locations of the contact. The location of the contact may be 
determined by plotting each hyperbola defined by each path length difference and 
selecting a point at which a large number of the hyperbolae intersect or nearly 
intersect. This point is likely to be the true location of the contact. 

[1018] Where n lm is unknown, the minimum number of series of hyperbolae 
required to determine the contact location is three and the likelihood of determining 
the correct location of the contact is increased by increasing the number of hyperbolae 
to be plotted. Multiple sensors may be used whereby a phase angle difference may be 
calculated for each pair of sensors thus generating multiple hyperbolae. In this 
embodiment, the minimum number of sensors is three. 

[1019] Alternatively, where ni m is unknown, the measured bending wave signal 
from each sensor may be divided into two or more discrete frequency bands whereby 
a phase angle difference may be calculated for each frequency band and for each pair 
of sensors. Although multiple phase angle differences may be calculated from a 
single pair of sensors, the phase angle differences at different frequencies are derived 
from the same path length difference. Thus the minimum number of sensors is three. 
The dividing of the frequency bands may be achieved by processing the bending wave 
signals by at least two band-pass filters having different pass-band frequencies. For 
example, using two band-pass filters having frequencies ooo+cosand g>o-cd& the phase 
angle differences A0 O , A0 b from two sensors may be defined as 

A0 a = k(co 0 + co 5 )Ax + 27m a 

A0 b = k{co Q - co s )Ax + 27m b 

where Ax is a single path-length difference defined by the contact and the position of 
the sensors. 

[1020| Therefore the values of n a and n b may be selected so that the measured 
phase angle differences infer similar values of the path-length difference. There may 
be only one combination of values (n a , n b ) for which this is possible. In this case the 
true value of the path-length difference may be determined. The correct combination 
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(n a , nb) may be determined as the combination of values that minimise the expression: 

A0 a -2m a A0 b -2m b 
k(co 0 +co s ) k(co 0 -a> s ) 

The path length difference may then be estimated as: 

^ = l( A0 a -2m a | A0 b -2m b ^ 
2{k{a> 0 +(a $ ) k(co 0 -co s )j 

[1021] Where this process is repeated with two pairs of sensors, two path length 
differences may be determined, which in turn may be used to determine the location 
of the contact. 

[1022] Alternatively, where ni m is unknown, an initial determination of the 
location of the contact may be made using the methods taught in WOO 1/48684 and 
PCT/GB2002/003073 (as summarized in Figure 1 1). Thereafter, if the condition 
Aco » 2k(a> 0 holds, the phase angle differences change by small increments over 

the timescale At = 2n I co . Accordingly, each value of n may be chosen to minimize 
the change in path length difference. 

[1023] Measured phase angle differences may contain random errors which may 
result in the selection of incorrect values of n. This error may be alleviated by 
evaluating the likelihood of successive sequences of n, for example by a state-space 
estimator such as the well known Kalman filter. The sequence having the maximum 
measure of likelihood is selected. 

[1024] The state-space estimator provides an estimate of the internal state of a 
system of which noisy measurements are made. A necessary input to the state-space 
estimator is a statistical description of the evolution of the system state. An example 
of such a state is the set of coordinates that describes the position and velocity of an 
object in contact with the member. It is widely known that the Kalman filter and other 
state-space estimators may provide a measure of likelihood that a sequence of 
observed, noisy measurements are consistent with the model of the system state. 

[1025] A state-space estimator may therefore be employed to take a sequence of a 
pair of path-length differences (say Ax\ 2 and Axr 34 ) taken at different times (say t\ 9 1 2 , 
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*3,...) s to estimate the system-state, i.e. the position and velocity of the contact, at 
those times. Moreover, the overall likelihood of those values of path-length difference 
being consistent with the model of the system may be evaluated. 

[1026] Where the sequence of path-length differences are obtained from a 
sequence of phase-angle differences and a set of integers {n = n(t\), n(t 2 ), wfo), . . ) 5 the 
measure of likelihood generated by the state- space estimator may be used to infer the 
likelihood that the correct values of n were chosen. It follows that a method for 
choosing the correct sequence of integers, w, is to find the sequence for which the 
state-space estimator gives the maximum measure of likelihood. 

[1027] As mentioned above, the state space estimator uses some statistical 
description of the evolution of the system state. A suitable model for the motion of the 
contact may be a simple random walk. Alternatively, the model may employ a 
detailed statistical description of how the user moves the stylus or finger. One 
example is a statistical description of how the user moves a pen while writing text or 
individual characters. 

[1028] The processor may further be adapted to include in the determination 
procedure any available information about where the contact can be expected. For 
example, if the member is an input device for a graphical user interface where the user 
is presented with a choice of 'buttons' to press, it may be useful to assume that any 
contact on the member occurs within the discrete areas corresponding to the buttons. 

[1029] Alternatively, a map of the probability at which a contact is likely to occur 
and which is based on the expected behaviour of the user may be used. The device 
may comprise a software application with a graphical user interface (GUI) which 
interacts with the operating system by means of an application program interface 
(API) in which the API is adapted to generate the probability map. The probability 
map may be based on the location, size, and frequency of use of objects presented by 
the graphical user interface. The probability map may also be based on information 
about the relative likelihood of the various GUI elements being activated. 

[1030] The following characteristics may apply to all embodiments of the 
invention. The device may comprise means for recording measured bending wave 
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signals from the or each sensor over time as the contact moves across the member. 
The information relating to the contact may be calculated in a central processor. The 
sensors may be mounted at or spaced from an edge of the member. The sensors may 
be in the form of sensing transducers which may convert bending wave vibration into 
an analogue input signal. 

[1031] The member may be in the form of a plate or panel. The member may be 
transparent or alternatively non-transparent, for example having a printed pattern. 
The member may have uniform thickness. Alternatively, the member may have a 
more complex shape, for example a curved surface and/or variable thickness. 

[1032] The device may be a purely passive sensor with the bending wave 
vibration and hence the measured bending wave signals being generated by an initial 
impact or by frictional movement of the contact. The contact may be in the form of a 
touch from a finger or from a stylus which may be in the form of a hand-held pen. 
The movement of a stylus on the member may generate a continuous signal which is 
affected by the location, pressure and speed of the stylus on the member. The stylus 
may have a flexible tip, e.g. of rubber, which generates bending waves in the member 
by applying a variable force thereto. The variable force may be provided by tip which 
alternatively adheres to or slips across a surface of the member. As the tip moves 
across of the member a tensile force may be created which at a certain threshold, 
causes any adhesion between the tip and the member to break, thus allowing the tip to 
slip across the surface. The bending waves may have frequency components in the 
ultrasonic region (>20 kHz). 

[1033] The member may also be an acoustic radiator and an emitting transducer 
may be mounted to the member to excite bending wave vibration in the member to 
generate an acoustic output. The frequency band of the audio signal of the transducer 
preferably differs from and does not overlap the frequency band of the measurements 
from the sensors. The audio signal may thus be filtered, for example, the audio band 
may be limited to frequencies below 20kHz, and the vibration measurements may be 
limited to frequencies above 20kHz. A sensor may have dual functionality and act as 
the emitting transducer. 
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[1034] The or each emitting transducer or sensor may be a bender transducer 
which is bonded directly to the member, for example a piezoelectric transducer. 
Alternatively, the or each emitting transducer or sensor may be an inertial transducer 
which is coupled to the member at a single point. The inertial transducer may be 
either electrodynamic or piezoelectric. 

[1035] A contact sensitive device according to the invention may be included in a 
mobile phone, a laptop or a personal data assistant. For example, the keypad 
conventionally fitted to a mobile phone may be replaced by a continuous moulding 
which is touch sensitive according to the present invention. In a laptop, the touchpad 
which functions as a mouse controller may be replaced by a continuous moulding 
which is a contact sensitive device according to the invention. Alternatively, the 
contact sensitive device may be a display screen, e.g. a liquid crystal display screen 
comprising liquid crystals which may be used to excite or sense bending waves. The 
display screen may present information relating to the contact. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[1036] The present invention may be better understood, and its numerous objects, 
features, and advantages made apparent to those skilled in the art by referencing the 
accompanying drawings in which: 

[1037] Figure 1 is a schematic plan view of a touch sensitive device according to 
one embodiment of the invention; 

[1038] Figure 2 is a schematic perspective view of the device of Figure 1; 

[1039] Figure 3 is a schematic side view of a one-dimensional beam; 

[1040] Figure 4a is a graph showing the amplitude of the reflection coefficient 
against frequency (Hz), the amplitude is unitless since it is a ratio; 

[1041] Figure 4b is a graph showing the phase (in radians) of the reflection 
coefficient against frequency (Hz); 

[1042] Figures 5a and 5b are schematic perspective views of alternative touch 
sensitive devices; 
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[1043] Figure 6 is a flowchart of a method of finding the location of a contact 
according to an embodiment of the invention; 

[1044] Figure 7a is a schematic block diagram of apparatus used for calculating 
phase angles; 

[1045] Figure 7b is a schematic block diagram of apparatus used with that of 
Figure 7a; 

[1046] Figures 8a to 8d are plan views of apparatus according to an embodiment 
of the invention showing the hyperbolae of path length differences; 

[1047] Figure 9 is a schematic block diagram of alternative apparatus used for 
calculating phase angles; 

[1048] Figure 10 is a flow chart showing an alternative method of calculating the 
location of the contact; 

[1049] Figure 1 1 is a flow chart showing a method of calculating the location of 
the contact using the dispersion corrected correlation function; 

[1050] Figure 1 la is a graph of dispersion corrected correlation function against 
time; 

[1051] Figure 12a is a schematic block diagram of a contact sensitive device 
which also operates as a loudspeaker, and 

[1052] Figure 12b is a method of separating audio signal and measured bending 
wave signal in the device of Figure 12a. 

[1053] The use of the same reference symbols in different drawings indicates 
similar or identical items. 

DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 

[1054] Figure 1 shows a contact sensitive device 10 comprising a transparent 
touch sensitive plate 12 mounted in front of a display device 14. The display device 
14 may be in the form of a television, a computer screen or other visual display 
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device. A stylus 18 in the form of a pen is used for writing text 20 or other matter on 
the touch sensitive plate 12. 

[1055] The transparent touch sensitive plate 12 is a member, e.g. an acoustic 
device, capable of supporting bending wave vibration. As shown in Figure 2, four 
sensors 16 for measuring bending wave vibration in the plate 12 are mounted on the 
underside thereof. The sensors 16 are in the form of piezoelectric vibration sensors 
and are mounted one at each corner of the plate 12. At least one of the sensors 16 may 
also act as an emitting transducer for exciting bending wave vibration in the plate. In 
this way, the device may act as a combined loudspeaker and contact sensitive device. 

[1056] In the following applications, U.S. Patent application number 09/746,405, 
filed 12/26/2000 entitled "Contact Sensitive Device" naming Nicholas P. R. Hill as an 
inventor, International Publication Number WO 01/48684 (International Application 
Number PCT/GB00/04851) and International Application PCT/GB2002/003073, filed 
July 3, 2002, which applications are incorporated herein by reference, contact 
sensitive devices and methods of using the same are described. The applications 
describe a device that includes a member capable of supporting bending wave 
vibration and a sensor mounted on the member for measuring bending wave vibration 
in the member and for transmitting a signal to a processor whereby information 
relating to a contact made on a surface of the member is calculated from the change in 
bending wave vibration in the member created by the contact. 

[1057] By bending wave vibration it is meant an excitation, for example by the 
contact, which imparts some out of plane displacement to the member. Many 
materials bend, some with pure bending with a perfect square root dispersion relation 
and some with a mixture of pure and shear bending. The dispersion relation describes 
the dependence of the in-plane velocity of the waves on the frequency of the waves. 

[1058] Bending waves provide advantages, such as increased robustness and 
reduced sensitivity to surface scratches, etc. However, bending waves are dispersive 
i.e. the bending wave velocity, and hence the "time of flight", is dependent on 
frequency. In general, an impulse contains a broad range of component frequencies 
and thus if the impulse travels a short distance, high frequency components will arrive 
first. U.S. Patent application number 09/746,405, International Publication Number 
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WO 01/48684 and International Application PCT/GB2002/003073, a correction to 
convert the measured bending wave signal to a propagation signal from a non- 
dispersive wave source may be applied so that techniques used in the fields of radar 
and sonar may be applied to detect the location of the contact. 

[1059] A mounting 22 made of foamed plastics is attached to the underside of and 
extends substantially around the periphery of the plate 12. The mounting 22 has 
adhesive surfaces whereby the member may be securely attached to any surface. The 
mechanical impedance of the mounting and plate are selected so as to minimise 
reflections of bending waves from the plate edges. 

[1060] The relationship between mechanical impedance of the mounting and the 
plate may be approximated by considering the one dimensional model shown in 
Figure 3. The model comprises a waveguide 34 in the form of a beam which 
terminates at an edge mounting 36 having a termination impedance. An incident wave 
38 travelling down the waveguide 34 is reflected by the mounting 36 to form a 
reflected wave 40. The incident and reflected waves are plane waves travelling in the 
direction perpendicular to the edge. Assuming the mounting 36 satisfies the 
following boundary conditions: 

(i) the termination impedance only couples into the lateral velocity, i.e. it does not 
provide any torque resistance; whereby the bending moment is equal to zero at the 
edge and 

(ii) the ratio of the lateral shear force and the velocity at the edge is equal to the 
terminal impedance; 

the reflection coefficient at the mounting is given by: 

Z T /Z B {co)+\ 

where Z T is the termination impedance of the mounting and Z B is the mechanical 
impedance of the end of the waveguide, given by 
_ / x Bk 3 (co), a 

Z *H= 2 ^ 0+0 

where £(co), is the wavevector which may be expressed in terms of the bending 



Final patent application 081-003.doc 
Client Reference: NXT 222 



- 12- 



Attorney Docket No.: 081-0003 

stiffness, B, and mass per unit area, jj,, of the panel, 




[1061] Thus, the reflection coefficient is determined by the ratio of the 
impedances at the end of the waveguide and the mounting. Furthermore, the 
impedance of the waveguide is proportional to the square root of frequency and is 
both real and reactive in equal weights (i.e. n/4 phase angle). Accordingly, the 
reflection coefficient is likely to be strongly frequency dependent. 

[1062] The reflection coefficient vanishes, i.e. bending wave energy is strongly 
absorbed in a frequency band around coo, if the following condition is satisfied: 
Z T =-iZ B (a) 0 ) 

[1063] Thus, the termination impedance of the mounting must have both real and 
imaginary components, or, equivalently, the mounting should be both resistive and 
compliant. 

[1064] The plate may be, for example, 1 mm thick polycarbonate sheet which has 
mass per unit area, /x=1.196 kg m" 2 and bending stiffness, B=0.38 N m. The 
equations above can be used to calculate the impedances of the plate and absorber 
required to strongly absorb bending wave energy around the chosen angular 
frequency G) 0 = 27r(900 Hz). 

[1065] The impedance, per unit width for a 1mm beam approximation of the plate 
is 

Zb(w 0 ) = (1+/) 33.8 Nsm" 2 . 

[1066] The properties of the absorber which provide the desired absorption are 
thus: Resistance per unit width, 

Re(Z r ) = Jm[Z B {co Q )] = 33 . 8 N s m" 2 . 
Stiffness per unit width, 

-ilm(Z r )co 0 =Re[Z B {co 0 )]cD 0 = 1.91 x 10 5 Nm' 2 . 
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[1067] The reflection coefficient is a unitless complex number. Figures 4a and 4b 
are graphs showing the amplitude and phase of the reflection coefficient R(co) varying 
with frequency. The amplitude of the reflection coefficient is zero and its phase is 
reversed for coq approximately equal to 900Hz. 

[1068] In Figures 5a and 5b, the plate 12 has uniform surface roughness in the 
form of a raised surface pattern 28,29. The stylus 18 is drawn across the surface along 
a path 30 and as it crosses a raised part or line of the pattern it generates bending 
waves 32 in the member. Thus contact from the stylus 18 provides a source of 
bending wave vibration in the member. In Figure 5a 5 the surface pattern 28 is a 
periodic pattern of raised crossed lines and in Figure 5b, the surface pattern 29 is a 
random relief pattern. 

[1069] In the embodiments of Figures 2, 5a and 5b, as the contact moves over the 
rough surface of the member, bending waves radiate isotropically in the member from 
the point of contact. The displacement of the member at a distance, x, from the point 
of contact is related to the displacement at the point of contact by a transfer function, 
H(co; x). At distances larger than the wavelength, X = 2n I k(a>) 9 the transfer function 
can be approximated as, 

ylk(a>)x 

where A is a constant and &(cd), is the wavevector defined previously. Although H(ca; 
x) strictly only applies to bending waves on an infinite plate, since the mounting 
strongly absorbs bending wave vibrations, the relationship is satisfied. The transfer 
function shows that where a source of bending waves emits a purely sinusoidal 
frequency with angular frequency, coo, the phase difference A0/2 between 
displacements at two locations which are at distances, x\ and x 2 , from the point of 
contact for the source is: 

exp(/A0 )2 ) = exp[i*(fi> 0 X*, - * 2 )] 
This implies the following relationship between the phase angle difference, the path 
length difference Ax = (x^ -x 2 ) and an integer ni 2 . 

A0 l2 =0 ] -6 2 = k(co 0 )Ax ]2 +273?i, 2 
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[1070] Figure 6 shows the steps in the method for using this equation to determine 
the contact location: 

(a) Measure a bending wave signal with each sensor to give measured bending wave 
signals W ( (t) and Wj(t), 

(b) Calculate the phase angles 0 g (t) and 0,(0 of the measured bending wave signals 
W.{t) and Wj(t), 

(c) Calculate the difference between the two phase angles 0.(f) and 0.(0 

(d) Calculate the location of the contact from 

k(a> o )Ax iJ =A0 iJ -2mi iJ 

[1071] Figure 7a shows a schematic block diagram of a device for calculating the 
phase angle 6j of a bending wave signal Wj(t) measured by one of the sensors. The 
signal Wj(t) is a random signal and is thus uncorrected over long time scales. The 
signal is first amplified by an amplifier 42 and then processed by an analogue band- 
pass filter 44 with a pass-band centered at coo and a bandwidth of Aco. 

[1072] A moving source of bending waves may demonstrate the Doppler effect, 
whereby a bending wave which has frequency coo and is emitted by a source moving 
at velocity v towards a point on a member arrives at that point with a different 
frequency defined by co 0 - k(co 0 )v . The maximum angular frequency shift between 

bending waves at two different points on the member is therefore 2t(® 0 )v nB , where 

v max is the maximum velocity of the moving source. If the angular frequency shift 
becomes larger than the width of the band pass filter, the phase difference equation 
above does not hold. Accordingly, the bandwidth Aco of the filter 44 is set to be 
greater than this maximum frequency shift and thus obeys the relationship: 
Aco » 2*(o> 0 )v max 

[1073] After processing by the filter 44, the resulting filtered signal Wj(t) is an 
amplitude and phase modulated carrier with frequency coo and is defined by: 

»^(f)=^(f>inK' + ^/0)J 
where Aj(t) and 0j(t) are the amplitude and phase of the signal. Both fluctuate over a 
timescale At determined by the bandwidth of the filter, namely At = 2n I Aco jh e 
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maximum frequency at which independent phase angle measurements may be taken 

from the output of the bandpass filter is . Since a touch sensor typically provides 

an updated measurement of the contact position every 10ms, the condition for the 
minimum frequency of positional measurement is At<10ms. 

[1074] The filtered signal W'j(t) is then passed simultaneously to two analogue 
phase detectors 46. Such detectors are well known in the art, for example, see p644 
of "The Art of Electronics" by Horowitz and Hill. Reference signals each having 
frequency coobut a phase difference of n/2 are also fed to the two phase detectors. 
The outputs of the phase detectors are passed through low-pass filters 48 each having 
frequency cut-offs of approximately Aca/2. The outputs of the low-pass filters are 
proportional to co%{0j) and s\n(0j) respectively. These outputs are then digitized by 
digitizers 50 and processed by processor 52 to give the phase angle Bp 

[1075] Figure 7b shows how the reference signals used in Figure 7a may be 
generated. A second bending wave signal Wi(t) is measured at a second sensor. The 
signal is fed through an amplifier 42 and analogue band-pass filter 44 to generate a 
filtered signal W\(t). The filtered signal W'tft) forms the reference signal which is fed 
directly to one phase detector 46. The filtered signal is also fed to the second phase 
detector 46 via a device which shifts its phase by n/2. The phase shifted signal is used 
as the reference signal to the second phase detector 46. 

[1076] Figures 8a to 8d show how the phase angle differences and hence the path 
length differences may be used to calculate the location of the contact. The equation 
in step (d) of Figure 6 defines a hyperbolic curve which can be overlaid on the plate 
12. Figure 8a shows the three hyperbolic curves 26 which are generated using three 
different values of ni m and the calculated phase angle difference for a pair of sensors 
16 mounted one on each end of the short sides of the plate 12. Similarly Figures 8b 
and 8c show the hyperbolic curves 26 which are generated by the phase angle 
difference and different values of ni m for two other pairs of sensors. Figure 8d shows 
all the hyperbolic curves created by the sensors. The contact location 24 is the point 
of intersection of three hyperbolic curves, one from each pair of sensors. From the 
contact location 24, the correct value of ni m may be inferred. 
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[1077] A method of inferring n is implemented using the embodiment shown in 
Figure 9. The bending wave signal W x (t) measured by each sensor is simultaneously 
processed by two band-pass filters 48,54. Two phase angles, one for each filter, are 
calculated, for example as described in Figure 7. The filters 48, 54 have slightly 
different pass-band frequencies whereby two phase angle differences, one for each 
pass-band frequency, are provided by each pair of sensors. 

[1078] The phase angle differences A0 O , A0b from the sensors may be defined as 

A0 a = k(co 0 + co 5 )Ax + 2rcn a 

A0 b = k(co 0 - co 5 )Ax + 2mi b 

where Ax is a single path-length difference defined by the contact and the position of 
the sensors. 

[1079] The correct combination (n a , n D ) may be determined as the combination of 
values that minimise the expression: 

A0 a -2jm a A0 b -27tn b 
k(a) 0 +co 3 ) k(co 0 -co s ) 

The path length difference may then be estimated as: 

A0 a -27m a | A0 b -2m b ^ 
k(a> 0 +co s ) kfa-ca^j 

[1080] Another pair of sensors may then be used to determine a second path 
length difference. Each path length difference defines a hyperbolic curve on the 
panel. The intersection point of these two hyperbolic curves is the location of the 
contact. 

[1081] Note that hyperbolae are defined by values of path length difference or Ax 
. In general, for a given phase-angle difference, several values of Ax are possible 
(corresponding to different values of n). The advantage of using two frequencies is 
that a single value of Ax can be obtained for each pair of sensors (using the method 
of minimizing the expression described above). The determination of the exact value 
of Ax , rather than a series of possible values, constrains the location of the contact to 
a single hyperbola, rather than a series of hyperbolae. The location can be determined 
exactly from the intersection of two hyperbolae, and hence from two pairs of sensors. 
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[1082] Figure 10 shows an alternative method for calculating the location of the 
contact from the equation above, namely 

i. Measure a pair of bending wave signals fV f (t) and W } (t) 9 one signal being 
measured by a sensor; 

ii. Calculate the dispersion corrected correlation function of the two signals using the 
method described in Figures 1 1 and 11a; 

iii. Calculate the initial position of the contact using the dispersion corrected 
correlation function, as described in Figures 1 1 and 1 la; 

iv. Remeasure bending wave signals W { {t) and W.(t); 

v. Calculate the phase angle of each signal - for example as described in Figures 7a 
and 7b; 

vi. Calculate the difference between the phase angles; 

vii. Select the value of ni m which minimizes the change in the path length difference; 

viii. Plot the hyperbola defined by 

k{w 0 )Ax iJ = AO-j -27m 0 

ix. Repeat steps (iv) to (viii), remeasuring the bending wave signals at regular 
intervals At, for example At =27c/Aco. 

[1083] At step (viii), a minimum of two hyperbolae from different pairs of sensors 
are required to determine the position of the contact. Therefore the entire is 
performed simultaneously for at least two pairs of sensors. 

[1084] Figure 1 1 shows a method of calculating the dispersion corrected 
correlation function to reveal the difference in path length between the contact 
location and the sensors. The method set out below summarizes the information in 
PCT/GB2002/003073. The method comprises the following steps: 

(a) Measure two bending wave signals W x (t) and W 2 (t) ; 

(b) Calculate the Fourier transform of W x (t) and W 2 (t) to arrive at W x (co) and 
W 2 (co) and hence the intermediate function W x (co) W 2 (co) ; where W 2 (co) is the 
complex conjugate Fourier transform, t represents time co is 2rcf where f is frequency. 

(c) Calculate a second intermediate function M(oo) which is a function of 
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(d) and (e) at the same time as performing steps (a) to (c), the frequency stretching 
operation f{co) = v{/u/B)^ 4co is calculated using the predetermined panel 
dispersion relation k = (/u/B)^ 4 4&> . 

(f) M(co) and f(co) = v(ju/B)^ 4 4& are combined to arrive at the dispersion corrected 
correlation function: 

G(t) = — (™M[f(6))]exp(iG)t)da> ; and 

(g) the dispersion corrected correlation function is plotted against time with a peak 
occurring at time t\ 2 as shown in Figure 11a; 

(h) Axn is calculated from t ]2 ; Axi 2 is the path-length difference between the path 
lengths xi and x 2 from the first and second sensors to the contact. 

(i) Ax 12 defines a hyperbolae which may be plotted as in Figure 7 to calculate the 
location of the contact. 

[1085] As with the method of Figure 10, a minimum of two hyperbolae are 
required to determine the location of the contact. Thus the ways of generating more 
hyperbolae discussed above apply equally to this method. 

[1086] The second intermediate function M(co) may simply be W x {co)W^ (co) 

which gives a standard dispersion corrected correlation function. Alternatively, M(go) 
may be selected from the following functions which all yield phase equivalent 
functions to the standard dispersion corrected correlation function: 



(a) M(ca) = 



(b) M(a>) = 



^\w,{co)w;(cd)\ 

(c) M(co) = ^(ty)^;(^(^(^)^;(^)|J where <p{x) is a real valued function 

(d) M(co) = W x {a>)W*{a>)\y{ai) where y/{co) is a real valued function 

[1087] Alternatively, M(o>) may be the function D(a>) which is the Fourier 
transformation of the correlation function D(t): 
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[1088] The steps are calculate D(t); calculate D(a>) and apply a frequency 
stretching operation to arrive at the dispersion corrected correlation function: 

G(t)= -I- ^D[f(co)]exp(ia)t)do) . 

[1089] Alternatively, at step (f) the following dispersion corrected correlation 
function may be calculated: G(t) = — PV, \f (co)}V 2 * \f{co)] <f> u \f{co)] exp(icot)dco 
where 



012 (*>) = 



where {^,, and {^*2,y(^)} are the Fourier transformation and complex conjugate 
Fourier transformation of two measured bending wave signals {^^(0} \fV 2 J (t)\ 
and {Ax 7 j is the path-length difference. 

[1090] A sensor may act as both the first and second sensor whereby the 

dispersion corrected correlation function is an autocorrelation function. The 
autocorrelation function may be calculated applying the same steps for the dispersion 
corrected correlation function using W x (t) = W 2 (t) . 

[1091] Figure 12a shows a contact sensitive device which also operates as a 
loudspeaker. Figure 12b shows a method for partitioning the audio signal and 
measured signal into two distinct frequency bands so that the contribution of the audio 
signal to the processed measured signal is suppressed. The device comprises a 
member 106 in which bending waves are generated by an emitting transducer or 
actuator 108 and the contact. The emitting transducer applies an audio signal to the 
member 106 to generate an acoustic output. Before being applied to the member, the 
audio signal is filtered by a low pass filter 1 12 which, as shown in Figure 12b, 
removes the audio signal above a threshold frequency f 0 . 
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[1092] As shown in Figure 12b, the contact generates a signal which has a power 
output which is substantially constant over a large frequency band. The signal from 
the contact and the audio signal sum to give a combined signal which is passed 
through a high pass filter 1 14 to remove the signal above the threshold frequency fo. 
The filtered signal is then passed to a digitizer 116 and onto a processor 118. 
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